Current-induced magnetization dynamics in disordered itinerant ferromagnets 
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Current-driven magnetization dynamics in ferromagnetic metals are studied in a self-consistent 
adiabatic local-density approximation in the presence of spin-conserving and spin-dephasing impu- 
rity scattering. Based on a quantum kinetic equation, we derive Gilbert damping and spin-transfer 
torques entering the Landau-Lifshitz equation to linear order in frequency and wave vector. Gilbert 
damping and a current- driven dissipative torque scale identically and compete, with the result that 
a steady current-driven domain- wall motion is insensitive to spin dephasing in the limit of weak fer- 
romagnetism. A uniform magnetization is found to be much more stable against spin torques in the 
itinerant than in the s-d model for ferromagnetism. A dynamic spin-transfer torque reminiscent of 
the spin pumping in multilayers is identified and shown to govern the current-induced domain-wall 
distortion. 



PACS numbers: 75.45.+j,72.25.Pn,72.15.Gd,72.25.Ba 

I. INTRODUCTION 

Metallic ferromagnets, notably the transition metals 
Fe, Co, and Ni, seem to be well understood, at least at 
temperatures sufficiently below criticality. Ground state 
properties such as cohesive energies, elastic constants, 1 : 
magnetic anisotropics in multilayers^ but also low- 
energy excitations that define Fermi surfaces* 3 - spin- wave 
dispersions, and Curie temperatures 4 are computed ac- 
curately and without adjustable parameters in the frame- 
work of local spin-density-functional theory (SDFT). 5 
Transport properties such as electric resistances due 
to random impurities are accessible to ah initio band- 
structure calculations as well^ However, important is- 
sues are still under discussion. Consensus has not been 
reached, e.g., on the nature and modeling of the Gilbert 
damping of the magnetization dynamics; 7 ^ the anoma- 
lous Hall effect^ and the current-induced magnetization 
dynamics. 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 The fundamental nature 
and technological importance of these effects make them 
attractive research topics. 

In this paper, we hope to contribute to a better 
understanding of the interaction of an electric current 
with a magnetization order parameter in dirty ferro- 
magnets, motivated in part by the sophistication with 
which the analogous systems of dirty superconductors 
has been mastered^ To this end, we proceed from 
time-dependent SDFT in an adiabatic local density ap- 
proximation (ALDA) and the Keldysh Green's function 



method in a quasiparticle approximation. We restrict 
ourselves to weak, diffusive ferromagnets where spin dy- 
namics take place near the Fermi surface, an approxima- 
tion that enables us to microscopically derive a simple 
quantum kinetic equation for the electronic spin distri- 
bution. The kinetic equation is used to derive a Landau- 
Lifshitz-Gilbert equation for the spatiotemporal magneti- 
zation that significantly differs from earlier phcnomeno- 
logical approaches based on the s-d model. We apply 
the general theory to the current-driven spin-wave ex- 
citation and domain-wall motion. After our paper was 
posted (cond-mat/0512715), Kohno et ali^S treated the 
same problem by diagrammatic perturbation theory. For 
weak ferromagnets, i.e., when the exchange potential is 
small compared to the Fermi energy, their results agree 
with ours. For strong ferromagnets, they report small 
corrections. 

The convincing evidence that transition-metal ground 
and weakly-excited states are well described by the mean- 
field Stoner model provided by local-SDFT can be ratio- 
nalized by the strong hybridization between the nearly 
free s-p bands and the localized d electrons^ It im- 
plies that the orbital angular momentum is completely 
quenched on time scales typical for the transport and 
magnetization dynamics. Both electric current and mag- 
netization are therefore carried by the same itinerant 
Bloch states. The alternative s-d model, in which only 
the localized d electrons are intrinsically magnetic and af- 
fect the delocalized s electrons via a local spin-dependent 



2 



exchange potential, is often used because it is amenable 
to sophisticated many-body treatments. On a mean-field 
level and with adjustable parameters, both models are 
completely equivalent for static properties. We find that 
the magnetization dynamics show drastic and experimen- 
tally testable differences that derive from the necessity 
of a self-consistent treatment of the exchange potential 
in itinerant ferromagnets that is not required in the s-d 
model. 

The paper is organized as follows: In Sec.lHl we discuss 
the model and the basic assumptions of the theory. In 
Sec. IIIII the quantum kinetic equation is derived in the 
real-time Green's function formalism, which is then used 
to obtain the magnetic equation of motion in Sec. IIVI 
The implications for the macroscopic dynamics are dis- 
cussed in Sec. [V] before the paper is briefly summarized 
in Sec. |Vl| 

II. MODEL 

In time-dependent SDFT) 21 i 22 i 23 the magnetic re- 
sponse is formally reduced to a one-body Hamiltonian 
in 2 x 2 Pauli spin space spanned by the unit matrix I 
and <x = (a x , a y ,a z ), the vector of one-half of the Pauli 
matrices: 

H = [H + U(r,t)+V[p](r,t)]l 

+ 7 ft«T-(H + H xc [p])(r,t)+H CT , (I) 

where Ti.o is the crystal Hamiltonian, U is the scalar dis- 
order potential including an external electric field, and 
V the spin-independent part of the exchange-correlation 
potential. We recognize on the right-hand side the Zee- 
man energy due to the sum of externally-applied and 
anisotropy magnetic fields H as well as an exchange- 
correlation contribution H xc , disregarding an exchange- 
correlation magnetic field coupled to the orbital motion. 
Here, 7 > is (minus) the gyromagnetic ratio and, H xc 
and V are functionals of the time-dependent spin-density 
matrix 

Po/s (r,t) = <¥j(r)* a (r)) t (2) 

that should be computed self-consistently from the 
Schrodinger equation corresponding to Ti. H. a is the spin- 
nondiagonal Hamiltonian accompanying magnetic and 
spin-orbit interaction potential disorder, thereby disre- 
garding the "intrinsic" spin-orbit interaction in the bulk 
band structure, apart from the crystal anisotropy con- 
tribution to H. Since we focus on low-energy magnetic 
fluctuations that are long range and transverse, we may 
restrict our attention to a single band with effective mass 
m e . Systematic improvements for realistic band struc- 
tures can be made from this starting point. We further- 
more adapt the ALDA form for the exchange-correlation 
field: 

7 /iH xc [p](r,i) « A xc m(r,f), (3) 



where m is the local magnetization direction with |m| = 
1 and A xc is the exchange splitting averaged over the 
unit cell. In terms of the spin density 

s(r) = HTr [&p(r)} , (4) 

m = — s/so, where Sq is the equilibrium value of |s|. For 
simplicity, the spin-independent random component of 
the potential U(r) is described as a zero-average, Gaus- 
sian white noise correlator: 

(U(r)U(r')) = £<5(r - r ') . (5) 

A characteristic scattering time r is defined by 



+ V \) T ' 

where v s is the spin-s density of states at the Fermi level. 
We consider two contributions to the spin-dephasing 
Hamiltonian spin-orbit scattering associated with 

the impurities and scattering at magnetic disorder that 
is modeled as a static random exchange field h(r) with 
white-noise correlator 

(h a (r)hp(r'))=t a 5 al3 8{r-r'). (7) 

It turns out both can be captured in terms of a prop- 
erly averaged, single parameter t„ for the characteristic 
transverse spin-dephasing time in the equation of motion 
for the magnetization. Derivation of the phenomenolog- 
ical T a for concrete microscopic models and dephasing 
mechanisms will be the topic of future correspondence. 

The ALDA is appropriate to describe corrections to the 
magnetization dynamics linear in d r that, although van- 
ishing for homogeneous systems^ are important in the 
presence of a current bias. The second-order correction 
(in homogeneous isotropic systems) is H' cx oc <r • <9 2 m, 
which contributes to the spin-wave stiffness [and can be 
taken into account via the effective field, see Eq. H29|) 
below]. Not much is known about the importance of 
nonadiabatic many-body corrections that in principle 
contribute to the magnetization damping. However, for 
slowly- varying perturbations of a homogeneous ferromag- 
net in time and space, the corrections to the ALDA are 
usually small^ Here we concentrate on dirty ferromag- 
nets in which the impurity (or phonon) scattering domi- 
nates quasiparticle scattering due to electron-electron in- 
teractions. 

In the next section, we derive the quantum kinetic 
equation for ferromagnetic dynamics by adiabatically 
turning on a uniform electric field until a steady state is 
established for a given current bias. The magnetization 
m is then perturbed with respect to a uniform ground 
state configuration mo = z. We then compute small de- 
viations of the spin density Ss = s + sqz, and replace 
s by — so m m the resulting equations of motion, com- 
pleting the self-consistency loop. A natural approach to 
carry out these steps is the Keldysh Green's function for- 
malism, which we briefly outline in the following. If the 
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reader is not interested in the technical details, we recom- 
mend jumping to Sec. II VI for the discussion of the result- 
ing equation of motion for the magnetization dynamics 
and Sec.[V]for the physical consequences for macroscopic 
dynamics. 



III. QUANTUM KINETIC EQUATION 

The Keldysh matrix Green's function can be repre- 
sented by the retarded G R {x,x'), advanced G A (x,x'), 
and Keldysh G K (x,x') components^ where x denotes 
position and time arguments. In the mixed (Wigner) 
representation (r, t; k, e), in which (r, t) are the center of 
mass coordinates, and using the gradient approximation 
(valid when hdt <C A xc and d r <C kp , a characteris- 
tic Fermi wave number), the Keldysh component of the 
Dyson equation reads in what is called the semiclassical 
approximation 



[G 1 , G K } p — {G K , G 1 ] p — 2i[G 1 , G K ] 



{"£ K , A} — {r, g k } 



(8) 



The left-hand side (l.h.s.) is the kinetic equation in the 
clean limit and the right-hand side (r.h.s.) is the collision 
integral. In the derivation of this equation, self-energy 
renormalization effects on the l.h.s and gradient correc- 
tions to the collision integral have been disregarded. This 
requires that A xc //i <C 1, where \i is the Fermi energy, 
although the corrections for large A xc appear to be very 
small, see below. £ is the self-energy due to disorder, 
which has three nontrivial components (i?, A, and K) 
along the Keldysh contour. Here, 



[B,C] p = d x B-d p C -d p B-8 x C 



(9) 



is the generalized Poisson bracket (where d x ■ d p = d r ■ 
9k — hdfd E ), [,] and {,} are matrix commutators and 
anticommutators, 



and 



A = i(G R - G A ) 



f = i{± R - t J 



(10) 



(11) 



G 1 is the inverse of the (retarded or advanced) Green's 
function in the absence of disorder: 

GoH'.SM) = [e-efc+e^(r,t)]i-Axc*-m(r,*), (12) 

where ip is the potential due to an applied electric field, 
and 



£k 



{hkf 

2m e 



/' 



(13) 



are the eigenvalues of Ho . We have disregarded the mag- 
netic field for the moment. In the self-consistent Born ap- 
proximation for scalar disorder scattering, the self-energy 
becomes 



E(r,t;k,e) =f / dk'G(r, t; k', e) 



(14) 



for each of the three components, where dk 1 = 
e? 3 k'/(27r) 3 . Self-energies for spin-dependent scattering 
channels can be calculated analogously. For A xc / [i -c 1 , 
we approximate the spectral function by Dirac delta func- 
tions at the two spin bands. Note that even though we are 
considering weak ferromagnets, the impurity concentra- 
tion is still considered dilute, so that H/t, H/t^ -c A xc , /i. 
By disregarding gradient terms of self-energies and the 
spectral function in the derivation of Eq. (JSJ , the Wigner 
representation transformed the collision integral into a 
local form. Gradient corrections disappear when the sys- 
tem is spatiotemporally homogeneous and/or we restrict 
our attention to weak ferromagnets, thereby discarding 
corrections of order O {ifi/r, A xc )/^). In spite of this re- 
striction, we believe that our formalism still captures the 
essential physics of the model (and therefore transition- 
metal ferromagnets) in a clear and coherent fashion. As- 
sessing the leading corrections to our treatment would 
require one to reconsider as well the simple ALDA mean- 
field treatment we are relying on. 

We concentrate now on the spin dynamics for small de- 
viations of the magnetization direction m = z + u from 
the z axis (u _L z) in the presence of a weak uniform elec- 
tric field E = —d r ip in the quasiparticle approximation 
for the Keldysh Green's function, 

G K (r, t; k, e) = -2ni ^ S(s - £ fcs )g ks (r, t) , (15) 

S 

where 



A 



(16) 



Two spin bands labeled by s =|, |= ± become separated 
when the disorder is weak. Note that in equilibrium, 



tanh 



£ks 



2k B T 



(17) 



where T is the temperature and ks the Boltzmann con- 
stant. The electric field applied to a rigidly- uniform fer- 
romagnet, u = 0, excites a nonequilibrium distribution 
3ks that is also diagonal in the spin indices. Interband 
spin-flip scattering vanishes upon momentum integra- 
tion, since a weak uniform electric field induces only a 
p-wave distribution. The transport in each spin band 
(obtained by integrating Eq. (|SJ) over energy e at fixed 
k near Sk s ) is thus described by the conventional Boltz- 
mann equation^ at T — > solved by the "drift" distri- 
bution 



Sgks 



he 



- + sa z E ■ VkS(eks 



(18) 



The distribution functions <?ks acquire off-diagonal 
components (describing transverse spins) in the presence 
of a finite u (so that out of equilibrium the spin subscript 
should not be taken literally). Equation |(HJ| leads to the 
linearized kinetic equation for the transverse component 
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gL = gks • <x (gks -L z): 

S^tgks + ^(vk • 9 r ) [g ks - A xc U(5(e fcs )] - A xc z x g ks 



sA xc z x usign(e fc;5 ) 



she 



(E • v k ) A xc z x ud(s ks 



•e(E. ak ) gks = 7 reE/*'^v-e fcs )[g k v- gks 
-(s - s')usign(e fcs )] + (v- s /v a - l)he (E • v k )u^(e fcs ) 



— g ks - su 

Tr 



he 

sign(e fcs ) + — ^E • v k <5(e fes ) 



(19) 



Quasiparticles propagate with group velocity v k = 
dk£k/h. On the l.h.s., an inhomogeneous exchange field 
is seen to cause electron acceleration and spin precession. 
The second term on the second line describes spin pre- 
cession of electrons accelerated by the electric field and 
the following term acceleration of the precessed electrons. 
On the r.h.s. we recognize elastic disorder scattering and 
transverse spin relaxation, the latter in terms of the spin- 
dephasing time t„. Energy-conserving mixing between 
the spin bands is allowed by disorder (in the presence 
of transverse fields), as reflected in the s' = —s part 
of the collision integral. We also took into account the 
contribution to the r.h.s. of Eq. ijSJ from anticommut- 
ing the current-induced drift Keldysh component with 
the spectral-function correction due to the magnetization 
deviation u: 



8 A = 2irfr 



u^]# - e ks ) ■ 



(20) 



IV. MAGNETIC EQUATION OF MOTION 

Integrating the kinetic equation l|19|l over momentum 
yields the equation of motion for the nonequilibrium spin 
density Ss = ~(h/4) J2 S J dkg Ws : 



r> r- A xc A xc 

Otds — z x os — z x uso 

h h 



^2 J dk(v k ■ <9 r )g ks 



Ss + USq 



(21) 



The integral on the r.h.s. is the divergence of the spin- 
current density, determined by the p-wave component of 
g ks , which can be found by a tedious (but straightfor- 
ward) manipulation of the kinetic equation. Confining 
our interest to spatially slowly varying phenomena re- 
sults in a major simplification: since d r already appears 
in Eq. I|21|l . we can disregard spatial derivatives in the p- 
wave component of g ks . We can now also include a static 
field H <C A xc along the z-axis by substituting primed 
quantities A xc = A XC +^H and u' = — (1— jH/ A' xc )5s/sq 
for the corresponding unprimed ones in the above expres- 
sions. The final result for the small- angle transverse spin 



dynamics is 
d t u = uj z x u — Pujqu+V 



1 — z x 



hd t 



(j-5 r )u, (22) 



disregarding the C(l/A xc ) terms inside the square brack- 
ets. Here j is the applied current density bias, 7H = luqZ, 







rA x 



(23) 



and V — (h/2e)P/so, where P = (cr-f — crj/^erj + o"x) i s 
the conductivity spin polarization, o s being the conduc- 
tivity for spin s along — m. For a Drude conductivity of 
parabolic bands, P = A xc /(ef'\+£fi)- We can transform 
the Bloch-like damping term in Eq. (122[l to the Gilbert 
form by multiplying the equation by 1 — (3zx from the 
left, which brings us to our central result: 



d t m = & ml 



LLG 



d t m\ 



where 



<9tm| LLG = —7m x H + /3m x d t m 



(24) 



(25) 



is the usual Landau-Lifshitz-Gilbcrt (LLG) with Gilbert 
damping 



0?LDA = P ■ 



d t m\ 



V 



1 - m x \/3 + 



hd t 



(26) 



(j-a r )m, (27) 



where, as before, we neglect the C7(l/A xc ) terms. 
«lda = h/r a A xc relates the collective magnetization 
damping to the single-electron spin relaxation that can 
be measured independently 8 and is consistent with ex- 
periments in permalloy films. 26 Equations H25JI and l|27|l 
hold for small deviations from a homogeneous equilib- 
rium state, but have the correct spin-rotationally in- 
variant form valid also for long-wavelength large-angle 
dynamics when the magnetic state is locally close to 
the equilibrium configuration (which requires a large ex- 
change splitting in comparison with other relevant energy 
scales). In particular, Eq. I|27() should correctly describe 
domain walls wider and spin-wave lengths longer than 
the magnetic coherence length hvp/A xc . For the same 
reason, the field H does not have to be nearly collincar 
with m. 

We can apply our method also to the mean- field s-d 
modelSI which leads to interesting differences. We repro- 
duced the phenomenologically derived Eq. (11) of Ref . IT?] 
(plus the dynamic term linear in <9 t ). The Gilbert damp- 
ing becomes reduced by the fraction 77 of the total spin 
angular momentum carried by the s electrons, while (3 is 
unmodified: 



(28) 
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assuming rj -C 1 (Ref . l27l) . We will see in the following 
that the ratio ft/a determines several interesting phys- 
ical quantities with /3/q;lda = 1 being a very special 
point. A sizable s-d character of the ferromagnetism al- 
ters this ratio, which could also be affected by a possible 
rf-magnetization damping in addition to the s-electron 
dephasing treated here. 

Soon after our paper was posted (cond-mat/0512715), 
a diagrammatic treatment of spin torques in static weakly 
disordered localized and itinerant ferromagnets has been 
reported by Kohno et al. 20 Their calculation is not re- 
stricted to weak ferromagnets [although it misses dy- 
namic current-driven torques such as the last term in 
Eq. J57J], and they find that in contrast to our result ft 
is not universally identical to a in the LDA approxima- 
tion. However, the ratio /3/q;ldAj m Ref. |20| expressed 
by the ratio between the density of states averaged over 
the two Fermi surfaces and the energy range spanned 
by A xc , is close to unity for almost all systems of in- 
terest. In particular, at low temperatures and in three 
dimensions, Kohno et al.'s expressions can be evaluated 
to be ft/ qilda ~ 1 + (1/48) (A xc /^) 2 (with the same 
correction for the ft/a s -d ratio). This quadratic devi- 
ation from unity is very small; even for A xc //i ^1/2 
it only amounts to about half a percent. The present 
quasiparticle treatment is not well suited to study ferro- 
magnets exhibiting an arbitrarily strong exchange split- 
ting (due to the increasing importance of gradient cor- 
rections to the semiclassical approximation with stronger 
exchange splittings). We are therefore hesitant to make 
predictions for half-metallic ferromagnets. We are, how- 
ever, confident that we capture the important physics of 
most experimental systems to date. For this reason, the 
present framework can also be used in studies of, e.g., 
relevant spin-dephasing mechanisms and microscopically 
derived scattering rates. The influence of realistic band 
structure effects, intrinsic spin-orbit, and Coulomb in- 
teraction, as well as corrections beyond the mean-field 
description might be more important than the gradient 
corrections to the ratio f3 /a. Furthermore, it is in general 
possible that other than impurity-related dephasing pro- 
cesses may contribute differently to a and (3, especially 
in the presence of strong anisotropies. 



V. CURRENT-DRIVEN DOMAIN- WALL 
MOTION AND BULK INSTABILITIES 

Let us proceed by discussing the influence of ft /a on 
the magnetization dynamics, and in particular the lim- 
iting case in which this ratio is unity. The dominant 
term r = V (j • d r ) m in Eq. H27|) is the conventional 
spin-transfer torque that, as far as the equation of mo- 
tion is concerned, can be absorbed into the magnetic free 
energyiii*i2ii£ The (dissipative) term proportional to ft 
acts like a magnetic field parallel to the direction of the 
magnetization gradient in the current direction. This 
term appears in our treatment by transforming Eq. I|22|) 



into the LLG form l|25|l . Zhang and Lii£ noted that al- 
though this "effective field" is much smaller than r when 
ft <C 1, it has a qualitative effect on the domain- wall 
motion. For example, in the absence of an external mag- 
netic field, a finite terminal velocity of a current-driven 
Neel wall is found for all currents only when the effec- 
tive field does not vanish. Judging from the importance 
of dynamic corrections to the spin torques in multilayer 
structures^ the dynamic contribution in Eq. Q27[) could 
be as significant since the typical frequencies of ferromag- 
netic dynamics are lu ~ r" 1 . 

In this section, we discuss several experimental conse- 
quences for j = jz, and a net effective field 

H= (Km z + H)z- K±m x x + AV 2 m. (29) 

Here, K is an easy axis and K± an easy-plane anisotropy 
constant, A is the exchange-stiffness, and H is the applied 
magnetic field. K, K±, A, H > 0. 

Let us first consider current-driven domain-wall mo- 
tion in the absence of applied field, H — 0. At the onset 
of the applied current density, a Neel wall along the z di- 
rection of width W = 'A/K with magnetization in the 
yz plane (pointing along z at z — > — oo and in the oppo- 
site direction at z — > oo) starts to movei£ with velocity 
(for not too large currents) 



(30) 



acquiring a terminal steady velocity for a constant cur- 
rent density given by 



v f 



ft 



Vj. 



(31) 



We find that the terminal velocity i|31|) is not influenced 
by the dynamic term on the r.h.s. of Eq. I|27|l . and we get 
Vf jvi = ft j q?lda = 1 for the self-consistent LDA model of 
itinerant ferromagnetism. The initial velocity l|3l)|l agrees 
with expectations based on angular-momentum conserva- 
tion, and, curiously, for our model, the terminal velocity 
is the same. According to Ref. |2(J, in three dimensions, 
the correction to ft/a^DA of order (A xc //x) 2 is positive, 
which means that Vf > Vi. Yamaguchi et al*& expressed 
the current-induced domain wall velocity 

v f = -(Vj 



(32) 



in terms of an "efficiency" ( of spin-current conversion 
into magnetization dynamics. Their experimental value 
C ~ 0.1 is much smaller than our result of £ = 1 in the 
absence of bulk or interface pinning (which, if smooth 
enough, could in principle be added to the effective field 
H) . For currents in excess of a threshold imposed by ex- 
trinsic pinning defects, Barnes and MaekawapS predicted 
C = 1 for an s-d model, in contrast to a nonuniversal 
mean- field result £ = ft/a s -d = 1/?/ of Ref. IT^ which we 
confirm here. 

Under the action of the current-induced spin torque, 
the shape of the moving domain wall distorts somewhat 
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with respect to the equilibrium configuration. The cor- 
responding domain-width change from the equilibrium 
value W to the steady-state value Wf was calculated in 
Ref. using the Walker's ansatz. After generalizing 
their method to include the effects of (3 as well as the 
dynamic term in the magnetic equation of motion (|27ll . 
we find 



w 



2-fA 



1 



(33) 



where the first (second) term on the r.h.s. describes 
the wall deformation due to the static (dynamic) part 
of Eq. ijTH) . Now, considering c*lda = (3, the first term 
vanishes and the wall slightly broadens, unlike the wall 
compression predicted for the s-d model with a finite 
damping a but setting (3 — (Ref. fl5^ . 

Finally, we discuss small-amplitude spin-wave solu- 
tions of Eqs. (JjSJl, |(2SJ, and 01 tne f° rm 



m(r, t) = z + Uo exp[i(q • r — cot)] 



(34) 



We are especially interested in solutions with Imw > 0, 
which describe exponentially growing spin-wave ampli- 
tude, signaling the onset of current-driven instabilities. 
We find that the critical current corresponding to Imw = 
is determined from 



a 



a A xc / V a A xc 



(35) 

where b' = P(q-j), H' = j{H+K+Aq 2 ), and K' = jK ± . 
For (3 — > 0, this reduces to 



\b'\ -> y/H'(K> + W 



(36) 



which can be thought of as the Doppler shift due to drift- 
ing spins necessary to overcome the natural spin-wave 
frequencyi 1:L i 12 i 14 i 15 Our result that oxda = (3 for weak 
ferromagnets, however, implies that a uniform magnetic 
state is stable against current-driven torques. In general, 
the critical current density j c determined from Eq. I|35|) 
can be significantly enhanced (depending on how close a 



and (3 are) with respect to the "Doppler-shift value" j c o 
calculated from Eq. I|36(l : 



3c 



JcO 



\l-(3/a\ ' 



(37) 



where small corrections proportional to (3 on the r.h.s. of 
Eq. have been disregarded. 



VI. SUMMARY 

In conclusion, we have used a quasiparticle approxima- 
tion, valid for weak ferromagnets, to derive an equation of 
motion for the magnetization dynamics of disordered fer- 
romagnets similar to the conventional LLG equation 1)25(1 
with Gilbert damping a and a current-induced contribu- 
tion (|27|l that is parametrized by a normalized single- 
electron spin-dephasing rate (3 = H/t^A^c. By virtue 
of the quasiparticle approximation, we obtain intuitively 
appealing kinetic equations that clearly reflect the phys- 
ical processes involved. 

Within a self-consistent picture based on the local den- 
sity approximation, we related the macroscopic damping 
in weak itinerant ferromagnets to the microscopic spin 
dephasing: chlda = f3, and pointed out striking implica- 
tions for current-driven macroscopic dynamics when the 
ratio j3/a is close to unity (which can also be expected for 
strong ferromagnets in the ALDA approximation). We 
furthermore noted remarkable differences in the dynam- 
ics of itinerant ferromagnets, supposedly well-described 
by the local-density approximation, and those with local- 
ized d or / electron magnetic moments. 
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